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Abstract 

A ratchet model for coupled Brownian motors, inspired by the motion of individual 
two-headed molecular motors on cytoskeletal filaments, is proposed. Such motors 
are modeled as two elastically coupled Brownian particles, each of which moves in a 
flashing ratchet potential. The ratchet potentials felt by the individual particles are 
anti-correlated, modeling the successive binding and unbinding of the two motor 
heads to the filament. We obtain, via Langevin simulations, steady-state currents 
as functions of noise strength, the equilibrium separation of the particles and the 
rate of switching between potential states. We observe an enhanced current due to 
coupling, noise induced stability and phase-locked behaviour in the deterministic 
regime. A qualitative understanding of these features is provided. 
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1 Introduction 



Motor proteins (kinesins, dyneins and myosins) are versatile biomolecules 
which shuttle cargo encapsulated in vesicles to different parts of the cell. Such 
transport occurs via the cytoskeleton, a cell-spanning polymeric network of 
microtubules, actin filaments and intermediate filaments. Motor proteins are 
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also important ingredients of the mechanisms of muscle contraction and cell 
division. The quest for physical principles that enable such tiny molecular 
machines to function efficiently in a highly Brownian regime is a subject of 
ongoing interest. [1] 

Molecular motors are typically dimeric[2]. Conventional kinesins, for 
example, consist of two identical proteins, each of which has a motor domain 
(the head) and a cargo binding domain. A stalk intervenes between the head 
and the cargo binding domain. The two stalks coil around each other to pro- 
duce the dimer. The two heads walk unidirectionally along a micro-tubule, a 
long and fairly stiff polymer comprised of an asymmetric repeating unit, in 8 
nm steps. Each such step is coupled to the hydrolysis of a molecule of ATP. 
Typically, such kinesins can walk several hundred steps before detaching from 
the micro-tubule track. 

Simple models for molecular motors idealize them as Brownian mo- 
tors[l,3]. A Brownian motor is a point particle with an overdamped equation 
of motion, which exhibits directed motion when subject to thermal and ather- 
mal noise [1,4]. The thermal noise is assumed to be Gaussian and satisfies 
a fluctuation-dissipation relation. For molecular motors, the athermal noise 
reflects the non-equilibrium driving of the motor through the irreversible hy- 
drolysis of ATP. 

The interaction of the motor with the molecular track is modeled in 
terms of an asymmetric periodic potential felt by the motor [1]. This potential 
switches between two or more states. The energy release upon ATP hydrolysis 
is believed to alter the coupling of the motor protein to the substrate. It is 
thus modeled by a switch in the potential state seen by the Brownian particle. 

The fact that the motor molecule is a two-headed object is ignored in 
these simple models (however, see Refs. [5,6]). Experimentally, two structural 
elements of the junction of the two heads appear to be crucial to rapid motility 
- the neck and the neck linker. The neck linker lies just beyond the catalytic 
core of the head, while the neck leads from the neck linker to the coiled stalk. 
The coupling of the two heads via the neck region ensures that an elastic in- 
teraction operates between the heads. The enhanced processivity and rate of 
movement in N. Crassa conventional kinesins has been attributed, in part, to 
the enhanced flexibility of its neck region[7,8]. 

Models of motor activity also indicate that the binding of one head as 
a consequence of ATP hydrolysis correlates closely with the loosening of the 
other head upon the subsequent release of the inorganic phosphate [9]. This 
suggests that a more physically realistic Brownian motor model of the work- 
ing of a two-headed molecular motor would involve (a) the incorporation of 
an elastic interaction between motor heads and (b) switching between poten- 
tial states which are fully anti-correlated between heads, so that one head is 
loosely bound when the other is tightly bound and vice versa. 

To develop these ideas further, we propose and study a simple model of 
a two-headed molecular motor interacting with a thermal bath at temperature 
T. Our model system consists of two elastically coupled overdamped Brownian 
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particles moving in periodic potentials Vi(x) and V 2 (x). These potentials al- 
ternate in time; we will assume this switching to be time periodic with period 
r. We take Vi(x) to be a periodic asymmetric potential and V 2 (x) to be a flat 
potential. The potentials felt by the two particles are out of phase i.e. when 
one of them moves under the influence of potential Vi(x), the other only sees 
the force due to the elastic linkage with the other particle. 

Our model has some similarities, and yet some important differences, 
with others proposed earlier. Unlike the models proposed by Klumpp et. al.[5], 
and by Ajdari[6], our model has no multiplicative noise terms, although the 
time-periodic variation of the potential does enter multiplicatively in our equa- 
tions of motion. Also, in our model, the potentials V\ or V 2 felt by one head 
are (anti-)correlated with the potentials felt by the other, not uncorrelated as 
in earlier work. 

We have studied our model numerically via Langevin simulations in 
the weak and strong coupling regimes. Our model exhibits a range of unusual 
behaviour, including an enhanced current due to coupling in some parameter 
range (this has been noted earlier in a different but related model [5]), drift 
in the deterministic regime (with phase-locked behaviour [11]), and noise in- 
duced stability [12]. Subsequent sections of this paper describe these features. 

In the next section we propose the system of coupled Langevin equa- 
tions to describe the model. In Sec. 3.1 we present results for our model in the 
deterministic regime, varying the coupling constant, the equilibrium separa- 
tion of the two particles and the frequency of potential fluctuation. In Sec. 3.2 
we discuss the behaviour of the current in the presence of thermal noise. Our 
conclusions are contained in Sec. 4. 



2 Model 

We consider two overdamped Brownian particles, coupled via an elastic 
interaction, moving in a two state flashing ratchet potential which switches 
between V\ (x) and periodically. The particles are in contact with a thermal 
bath at temperature T. Let x±(t) and X2(t) be the coordinates of the two par- 
ticles. They are coupled through a spring of spring constant k and equilibrium 
length a. The Langevin equations which governs our system are 



where £j(t) are Gaussian random forces with zero mean and a correlation 
given by < (f) >= 2£>(%5(t-t / )- Here, D = 4s2l. The ^'s are periodic 

functions with period r, given by Z\{t) = l,z 2 (t) = for < t < ^;zi(t) = 
0, ^(O = 1 for I < t < r. The asymmetric periodic potential with unit period 
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2i(0^'(xi) + k((xi - x 2 ) -a)+ £i(0 
z 2 {t)V'(x 2 ) - k(( Xl - x 2 ) -a)+ 6(0 



(1) 
(2) 
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is V(x) = — |^(sin(27r:r) + 0.25 sin(47rx)). We work with dimensionless units 
wherever possible. 

We solve these coupled equations numerically by Huen's method [13] 
and calculate the current in the asymptotic regime. The current j is given by 
3 — ( ^zf~ )) where < .. > denotes the ensemble average and t the asymp- 
totic time. We take to = 50r, t = 4000r and average over 3000 ensembles. In 
the limiting case of rigid rod like coupling (i.e., x\ — X2 = a), we solve the 
associated Fokker-Planck equation 9P q*'^ = — 9J q^ [14] 

^-|(V(.) + ^|)P(^), (3) 



for nr <t< (n+ l/2)r, 



dP(x.t) d , /— d . „. 



(4) 



for (n + l/2)r < t < (n + l)r. Here x is the position of one of the particle, 
the position of the other particle being fixed at x + a and P(x, t) is the proba- 
bility density. The average current j = lim^oo - j/ +r dt Jq 1 J(x,t)dx. We have 
obtained j by numerically solving the Fokker-Planck equation. 



3 Results and Discussions 



3.1 Deterministic Current (D = 0) 



A single particle moving in a flashing potential such as the one de- 
scribed above, exhibits no directional motion in the absence of thermal fluctu- 
ations. However, such motion can be induced in a system of elastically coupled 
particles using the compressibility of the spring [6,5]. If the equilibrium sep- 
aration of the two particles is smaller than the length of the smaller arm of 
the potential V(x) [L m in = 0.38) or larger than the long er arm {L m ax — 0.62), 
then the phase space of the system is a closed orbit. No directed current is 
possible for any value of r or k. The window L max — L m i n , for which a directed 
current is obtained depends on the value of r (Fig.(lA)). For smaller values 
of r, the length of the window for nonzero current is small, as the particle has 
less time to traverse to the basin of attraction of the next well. In Fig. (IB) 
we plot the current j as function of k for a = 0.5, for different values of r. 
Observe that the current exhibits a phase locked behaviour, with step size 
becoming larger for smaller r as observed previously in many different sys- 
tems. Such phase locking arises in deterministic systems due to the interplay 
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Fig. 1. [A] Current j vs equilibrium distance a at k = 0.8; [B] current 
j vs coupling constant k for a = 0.5. The inset shows j vs r for same parameter 
value. 

between nonlinearity and external drive. 

For large k, the current saturates to a constant value, which depends 
on r. At these values of k, the particles traverse exactly one spatial period in 
one time period; since our system is overdamped, no multiple hops can occur. 
The inset shows the variation of j with r for k = 3.0 and k = oo (rigid rod 
case). The value of r for which current is maximum decreases with increasing 
k. 



3.2 Effect of Noise 



Many of the attributes of the deterministic limit, such as the phase 
locked behaviour of Fig. (IB), are sensitive to noise. In the presence of ther- 
mal fluctuations the noise-induced current exhibits many interesting features. 
There are three distinct regimes in the parameter space where qualitatively 
different features are obtained: (i) large r and L min < a < L max , (ii) large r 
and a < L min and (iii) the small r or nonadiabatic regime, which we do not 
discuss here. 

We first study case (i), examining the effect of noise in a parameter 
range where there is current in the deterministic limit. In the large r, a = 0.5 
regime the current j increases initially with k for all values of D as shown 
in fig.(2A). On further increasing k (not shown) the current saturates after 
a minimum. The inset shows the variation of j with D for different values 
of k, keeping the other parameters the same. The case k = is the single 
particle limit. The location of the peak in the j — D curve is determined by 
the time taken to diffuse to the basin of attraction of the next well (j max at 
D l? min jT = 0.16/5 = 0.032). For non-zero coupling, when one particle 
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Fig. 2. [A] Current j vs k for a = 0.5. The inset shows the plot of j vs D for 
k = 0.0, k = 0.3, k = 0.5 and k = 0.67. ( From top to bottom curve.) [B] The peak 
in the j vs D ( corresponding to minimum current) emphasizes the noise induced 
stability phenomena. Here a = 0.5. 



slides down the potential the other is dragged along, thereby enhancing the 
current. However, on further increasing the noise strength, the potential land- 
scape becomes irrelevant, random hopping in both directions dominates and 
the absolute value of the current tends to zero, accounting for the minimum 
seen. At k = oo, the current decreases monotonically with increasing D. 

Though the increase of current with D for low coupling strength or 
monotonous decrease at higher k is the generic feature in this regime of op- 
eration, for certain values of k (values for which there is phase locking in the 
deterministic regime, see Fig. (IB)) the current shows highly nonintuitive be- 
haviour. As shown in Fig. (2B), the absolute current initially decreases with 
increasing D, reaches a minima and then again increases with D. This feature 
is reminiscent of the noise induced stability (NIS) of unstable states [12], seen 
previously in different models, i.e. noise localizes the particle in its basin of 
attraction. This leads to a decrease in the value of the current with noise. In 
the high temperature regime, the current should go to zero, explaining, qual- 
itatively, the two extremas seen in the current as a function of noise. 

We now investigate the effect of noise in the parameter range where 
there is no current in the deterministic limit, our case (ii). We work with 
a = 0.2. Fig. (3) shows j vs D for various values of k. Unlike the previous case 
where the temperature at which the current peaks decreased with increasing 
k, here this temperature first decreases and then increases with increasing k. 

The current here arises as a consequence of both the "flashing ratchet 
mechanism" and the cooperative "pulling effect". For low coupling, the con- 
tribution of coupling at a temperature at which the current arising out of 
the "flashing" mechanism is maximum, is negligible. Hence the net current 
peaks at a lower value of T compared to the k = case. For higher k, where 
the current due to coupling dominates, a different scenario obtains. As the 
stiffness of the spring is larger, a larger fluctuation is required for the parti- 
cle to transit between wells. The larger the k, the stronger the amplitude of 
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Fig. 3. j vs D for different values of k, a = 0.2. The inset shows the variation of 
j vs k in log scale, for D = 0.01, 0.04, 0.07and 0.1 from top to bottom. 

the fluctuation required. If, however, the temperature is too high, significant 
motion is induced in the opposite direction and the current reaches a maxi- 
mum. Depending on the parameter regime, the coupling constant may or may 
not enhance the current compared to the single particle current. This can be 
clearly seen from Fig. 3 and the inset. This behaviour again differs from the 
case where there is a finite current in the deterministic limit. 

In the nonadiabatic regime obtained for low values of r, even for ex- 
tremely small values of k, currents may either enhance or decrease relative to 
the single particle current. This differs from our earlier observation in the adi- 
abatic limit where an enhanced current is obtained for arbitrary small values 
of k, for any D. These results will be reported elsewhere. 



4 Conclusions 

Our simple model of a two-headed motor comprised of two coupled 
Brownian particles in a flashing potential, exhibits unidirectional current. The 
coupling between particles induces an additional mechanism by which a cur- 
rent may be obtained (via a pulling effect), which is absent in the single particle 
case [5]. The behaviour of the current depends sensitively on system parame- 
ters and exhibits several novel features, such as an enhanced current due to 
coupling compared to the single particle case, in a regime which we identify. 
Outside this regime, coupling may or may not lead to a current enhancement. 

As a function of noise strength, the current can exhibit noise induced 
stability in a restricted parameter range. Outside this range, noise initially 
enhances current and the absolute value of current shows a single maxima. 
In the deterministic case, we have shown the existence of a window for the 
two particle separation distance within which there is a deterministic current. 
This current shows a phase locked behaviour with coupling constant. Unlike 
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previous studies on related models, we have not observed current reversals as 
a function of system parameters. 

What relevance, if any, might these results have to real motors? While 
generically, weakening the coupling of the two heads (particles) appears to 
reduce the current, we do see specific regimes in which such a reduction leads 
to an enhancement in the current. As discussed in the Introduction, the in- 
creased motility of some Kinesins has been attributed to a larger flexibility 
of the neck linker. However, whether such simplified models can ever reflect 
the underlying behaviour of vastly more complex biological systems is unclear. 
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